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Abstract. The presence of ultra-high energy cosmic rays in the universe in the pre-ionization epoch enhances 
significantly H2 kinetics in virialized halos. It results in more efficient radiative cooling and decreases the lower 
mass limit of the first star-forming systems. Consequently, the fraction of baryons contained in the first luminous 
objects and their contribution to the reionization of the universe can significantly increase in comparison with the 
standard scenario. 
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1. Introduction 

In hierarchical cold dark matter (CDM) scenarios the 
formation of the first stellar objects involves dark mat- 
ter halos of total (dark and baryonic matter) masses 
Mh ^ 2 X 10*5 Mq at redshifts z 20 (Tegmark et. al. 
1997, hereafter T97), or lower {Mt ~ 2 x lO^M©) if the 
effects of violent relaxation of dark halos are taken into 
account (Vasiliev & Shchekinov 2003). This process, and 
the corresponding halo masses where the first stellar ob- 
jects may form, is determined by the cooling of baryons 
mainly through radiation of thermal energy in rotational 
lines of H2 molecules. In primordial gas, H2 molecules can 
form only in reactions involving H~ and ions, and 
therefore the H2 yielded strongly depends on fractional 
ionization x = Ue/n (Galli & Palla 1998). In the stan- 
dard cosmological recombination, H2 can reach not more 
than /2 = n(H2)/n = 10^"^ in already virialized halos and 
it determines a rather high lower mass limit Mh (T97). 
However, in cosmologies with decaying dark matter par- 
ticles, the recombination dynamics itself changes, as was 
first suggested by Sciama (1982, 1990), Adams, Sarkar & 
Sciama (1998). As a result, H2 chemistry can be signifi- 
cantly altered. 

Ultra- high energy cosmic rays (UHECR), if they 
formed in the so-called top-to-bottom scenario from 
decaying superheavy dark matter (SHDM) particles 
with masses Mx ^ 10"'^^ GeV (Berezinsky, Kachelriefi 
& Vilenkin, 1997, Kuzmin & Rubakov, 1998, Birkel & 
Sarkar, 1998), may also have a strong impact on cosmo- 
logical recombination. A detailed description of the ion- 
ization history of the primeval plasma in the presence of 
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such UHECR has been reported recently by Doroshkevich 
& Naselsky (2002, hereafter DN), Doroshkevich et al. 
(2003). These authors conclude that the UHECR con- 
verted through electromangetic cascades into UV photons 
increase the hydrogen ionization at z ~ 10 — 50 by a fac- 
tor of 5 to 10. Correspondingly, the fraction of molecular 
hydrogen in virializing halos can increase by factor of 10- 
20 depending on details of their dynamics. In this paper 
we discuss briefiy the effects of such enhanced kinetics of 
H2 on minimal masses and evolution of early luminous 
objects. 

Throughout this paper we assume a ACDM 
cosmology with parameters {flo,flA,^m,^b,h) — 
(1.0, 0.71, 0.29, 0.047, 0.72) as inferred from the 
Wilkinson Microwave Anisotropy Probe (WMAP), and a 
deuterium abundance of 2.62 x 10^^ (Spergel et al. 2003). 

In Section 2 we estimate minimal masses of early lu- 
minous objects. In Section 3 we analyze qualitatively how 
additional ionizing photons change thermal dynamics of 
individual contracting baryon condensations. In Section 4 
we discuss possible cosmological consequences. Section 5 
summarizes our conclusions. 



2. Minimum masses of early luminous objects 

Formation of early luminous objects (stars, quasars or stel- 
lar systems) in the universe and their masses are deter- 
mined by the ability of primordial gas to cool radiatively 
and to form sclf-gravitating contracting agglomerations of 
baryons. In the hierarchical CDM scenario, contraction 
of baryons is catalyzed by growing fluctuations of gravi- 
tational potential from virialized dark matter halos. The 
exact description requires a full 3D multi-fluid numerical 
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approach, and although such a description explains the 
whole variety of dynamical processes involved, it remains 
still highly time consuming (Abel, Bryan & Norman 2000, 
Bromm, Coppi & Larson 2001, Abel, Bryan & Norman 
2002). For our purpose, to illustrate possible consequenses 
of the presence of extra ionizing photons generated by 
UHECR on the process of the initial star formation, one 
can follow a simplified approach suggested by T97. This 
approach is based on the solution of an ordinary differ- 
ential equation for thermal energy of baryons confined by 
the gravitational potential of a virialized halo. An essen- 
tial assumption of T97 is in their criterion for baryons 
to cool: they assumed that the baryons in a halo with 
mass M will cool progressively, if the baryon temperature 
decreases sufficiently: T{riz) < rjT{z) faster than in one 
comoving Hubble time tniz), where r] = 0.75. 

Vasiliev & Shchekinov (2003) have used a different cri- 
terion which requires stronger conditions for baryons to 
cool and form a self-gravitating configuration: they as- 
sumed that it occurs if in the process of baryon contrac- 
tion in the gravitational field of a dark halo, the baryons 
reach a mean density equal to the mean density of dark 
matter within the sphere occupied by baryons. Obviously, 
it results in an increase of the minimal baryon mass able to 
cool. However, to qualitatively understand possible effects 
from UHECR on star formation, here we present calcula- 
tions of the minimal mass only within the T97 criterion. 

Consider first a single virialized halo of mass Mh = 
10^ Mq. The equations describing ionization fraction x = 
n(H+)/n, molecular fraction / = n(H2)/n and thermal 
evolution of baryons read as (T97) 



X = —kinx^, 

f = kmn{l -X- 2f)x, 

T= ^-T + Ac-Ahi-Ah,, 
3 n 



(1) 
(2) 



(3) 



where Ac, Ahi, Ah2, Ahd are the Compton, HI, H2 and 
HD cooling rates, respectively, ki and km are the reduced 
reaction rates (see for details T97). A reduced HD chem- 
istry has been also added to account for effects from HD 
cooling at temperatures lower than 400 K. We have run 
the system for a grid of halo masses Mh with initial con- 
ditions corresponding to the fractions of species and the 
temperature at the turnaround redshift; the evolution be- 
fore the turnaround point has been found from the T97 
model. 

The rate of production of ionizing photons by UHECR 

can be written as dn^/dt = e{z)H{z)n{z) where is the 
number density of ionizing photons, H[z) is the Hubble 
parameter and n{z) the gas number density (see DN). 
Here the efficiency e{z) is defined as (DN): 



2.5 X 10"^ 
l + z 



(4) 



where a is the power index of the spectrum of injected 
photons for a decay of a single SHDM particle, Miq = 
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Fig. 1. The ionization (solid), molecular (dash) and HD 
(dot) fractions, and the thermal evolution (dot-dash) of 
baryons in a virialized halo with mass Mh = 1O^M0. The 
upper panel depicts standard case without additional Lyc 
photons, while the lower panel represents a model with 
extra UV photons with parameters corresponding to the 
model e = 0.3/(1 + z) in Doroshkevich & Naselsky (2002). 
The HD fraction is multiplied by 10^. The number density 
n is divided by 10^. 

Mx/lO^'^ GeV, Mx is the mass of the SHDM particles. 
For the measurable distortions of the power spectrum of 
the CMB anisotropy Qtot - lO^Mie""'^ (DN). In the cal- 
culations we assume for the efficiency 0.3/(1 -I- z) as a 
typical value. 

In Fig.l the thermochemical evolution of baryons in a 
virialized halo with mass Mh = lO^M© for the standard 
recombination scenario (upper panel) and that with an 
enhanced fractional ionization (lower panel) is shown. The 
small additional electron fraction produced by UHECR 
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photons results in more efficient H2 molecule formation, 
such that inside a collapsing low-mass halo the molecular 
hydrogen fraction reaches a value of about (0.6 — 4) x 
10~^ depending on e. As a concequence the temperature 
drops faster, and this enhances formation of HD molecules, 
and practically all deuterium is bound into HD molecules 
when the temperature falls below 200 K. Thus, even for 
a relatively low production rate of ionizing photons e the 
molecular cooling rate increases significantly. 

In principle, additional UV photons from UHECR can 
heat baryons in the halo, however in most cases their im- 
pact on thermal evolution inside a formed halo seems to 
be unimportant: the corresponding photoionization heat- 
ing is Hjjv ~ 2 X lO^^^e n(l-l-z)^-^ erg cm^^ s"^ [n is the 
baryon density) and can be cooled by molecular hydrogen 
with n(H2)/n - 3 x lO"'' at T - 700 K, i.e. smaller than 
virial temperature for halos with Mh > 3 x IQ^Mq. This 
means that when a halo of Mh > 3 x lO^Af0 reaches its 
virial state with > 700 K, radiative energy losses by 
H2 molecules exceed UV heating. Note also that the rate 
destruction of ions by extra photons is only about 
10~^^s~^, which is much smaller than the rate of conver- 
tion to H2 lO^i^g-i 

In such conditions self-gravitating baryon configura- 
tions form more easily and have the lower masses. Fig. 2 
shows the z-dependence of lower mass limit Mhi{z) of ha- 
los able to form cooling self-gravitating baryon condensa- 
tions. The thick solid line depicts the masses correspond- 
ing to 3cr peaks in a ACDM model. Only masses in the in- 
terval Mhi{z) < Mh < M3cr{z) are supposed to form early 
luminous objects. In general, additional ionizing photons 
from UHECR can stimulate formation of early objects in 
the universe with masses of an order of magnitude lower 
than in conditions with the standard cosmological recom- 
bination. 

3. Evolution of individual halos 

Besides determining the minimal mass of the first baryon 
objects, the question of how individual halos of large 
masses evolve is also of interest for a qualitative under- 
standing of the effects of UHECR on primordial star for- 
mation. Of course, the full description can be reached only 
in multi-dimensional numerical simulations, however, to 
understand qualitatively what one can expect from the 
presence of UHECR in the early universe for primeval star 
formation, a simplified consideration may be helpful. 

Within this approach we will describe baryons in a 
virialized halo by the following set of equations widely 
used in simplified calculations (see e.g., Tegmark, Silk & 
Evrard 1993) 



R = u, 



u = 



T - 



finipR 3 
2kT 



T^G{pd + p)R, 



(5) 
(6) 

(7) 




1+Z 



Fig. 2. Minimum mass able to cool and collapse. The 
dashed line depicts the lower mass limit estimates by T97 
in standard cosmological recombination. The four thin 
solid lines represent the lower mass for the rate of ionizing 
photon production by UHECR, e (0.1, 0.3, 1, 3)/(l -h z) 
(from top to bottom). The thick solid line corresponds to 
3ct peaks in a ACDM model. 



where all notations have their usual senses, A^ are as in Q. 
These equations were solved with a complementary set of 
destruction for ionization and chemical kinetics account- 
ing for additional ionization of H and photodestruction of 
H^ ions by UHECR photons, and a simplified treatement 
of effects of optical thickness on the corresponding cooling 
function by introducing an exponential factor exp(— rpja), 
where is the average optical depth in H2 lines. We cal- 
culated the evolution of baryons for a halo with a given 
mass from its virialized state to the asymptotic state when 
concentrations of H2 and HD reach their limiting values, 
and temperature settles onto a regime determined by adia- 
batic heating, radiative cooling and energy exchange with 
cosmic background radiation. 

The results are shown in Fig. 3, where temperature, 
baryon density and Jeans mass within a baryon conden- 
sation are depicted, respectively for three halo masses: 
Mh = 10^, 10^ and IO^Mq. 

The difference between the temperatures is readily 
seen: in the case e = (standard recombination), after 
an initial period of cooling, temperature starts growing 
for the halo mass Alh = IO^Mq because strong adiabatic 
heating cannot be balanced by H2 with H2/H'^ 3 x 10^**, 
and only in subsequent contraction when H2/H increases 
fast and reaches 5 x 10~^ temperature decreases again. 
It happens only at z ~ 13.5. In the case of nonzero in- 
put from UHECR, due to enhanced production H2 abun- 
dance reaches sufficiently high values ^ lO^'^ already at 
the initial stages, so that a powerful H2-cooling sets in 
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Fig. 3. Temperature (upper panel), density (middle 
panel) and Jeans mass (lower panel) evolution of baryons 
in halos with masses 10^, 10^ and 10^ M© virialized at 
z = 20 for the standard scenario, e = 0, and for cosmolog- 
ical recombination dominaned by UHECR photons, e = 
0.1, 3: solid lines correspond to halo mass = W^Mq, 
dashed lines - Mh = IO'^Mq, dotted lines - Mu = 10^ Mq. 
On panel (a) the temperature of CMB is shown by the 
thick solid line. 



at the very initial contraction phase, and the tempera- 
ture decreases monotonously down to ~ 50 K (the CMB 
temperature) in a much shorter time - the corresponding 
redshift z ^ 16.5 for e > 0.3. At the latest phases of the 
contraction (T ^ 200 K), the thermal regime is controlled 
by cooling from HD molecules whose abundance reaches 
the abundance of deuterium in models with e ^ 0.3. Halos 
with masses = lO^M© show a similar behavior de- 
pending on e (as shown by solid lines in Figs. 3), although 
the difference between the models with e = and e > is 
less pronounced here. Thus, the effects of additional ion- 
ization are more important for the halos at the low end of 
masses. 

Generally, the Jeans mass within the baryonic core de- 
creases faster when the ionization efficiency from UHECR 
photons increases. At later stages it reaches asymptoti- 
cally the value Mj ~ 100(1 + z)^/^ n-^^M© while bary- 
onic cores or fragments remain optically thin for HD lines, 
here the dependence on z stems from asymptotical value of 
baryon temperature Tb ~ To(l -|- 2;). When the Jeans mass 
becomes opaque in HD lines, its value is 1-2 Mq. For the 
standard scenario (e = 0), the Jeans mass at later stages 
of contraction is nealy coincident (within factor of 2) with 
the mass of a collapsing core Mc ~ 100 Mq reached in 
numerical simulations by Abel, Bryan & Norman (2003). 
Note that in their simulation the temperature of gas does 
not decline to the CMB value - the typical temperature in 
the core is ~ 300 K. Our results differ from the results of 
a full 3D simulation apparently due to contribution from 
HD molecules. 

4. Discussion 

An obvious consequence of the enhaced H2 production dis- 
cussed here is closely related to the post-dark age evolu- 
tion of the universe: its reionization and enrichment with 
metals by the primary stellar nucleosynthesis. The pres- 
ence of an additional ionization from UHECR photons de- 
creases the lower limit for the first baryon condensations 
and shifts their formation to earlier epochs: with simple 
T97 criterion for formation of baryon objects, additional 
UV photons decrease the minimal mass from ~ 10^ Mq 
to ^ 10^ Mq, and correspondingly increase the redshift 
from z ~ 20 to z ~ 25. In a hierarchical scenario (e.g. 
Lacey & Cole, 1993) the number of 10® Mq halos at red- 
shift z « 20 is an order of magnitude smaller than of 
IO^Mq. This means that the production of Ly-c photons 
by the first stars may have been significantly greater. If 
the efficiency of star formation in smaller halos is equal 
to that in more massive halos, the production of photons 
increases by a factor of 2, which would be sufficient to 
reionize the universe at levels established by WMAP (see 
discussion in Cen 2003). However, if the efficiency of star 
formation depends on halo mass as /* ~ M", than with 
a = 0.8 as suggested for instance by Kasuya et al (2003) 
the amount of Ly-c photons increases only by ~ 20 %. 

The more numerous first baryon objects expected in a 
model with the presence of UHECR photons may change 
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significantly also the metal enrichment history of the inter- 
galactic medium (IGM), as observed in Ly-a forest clouds. 
One of the key problems, related to the initial pollution 
of the IGM by metals is connected not only with their 
total amount observed in Ly-a absorbers, but with their 
widespread presence at high redshifts, and particularly in 
underdense regions far from possible galactic sources (see, 
the recent discussion in Schaye et al 2003) . A possible so- 
lution of this problem can be related to the highly overlap- 
ping metal - enriched outflows produced by the early star- 
forming halos, as suggested by Nath & Trentham (1997), 
Madau et al (2001), Scannapieco et al (2002). In particular 
Scannapieco et al (2002) predict that halos of intermedi- 
ate masses of 3 x 10®^^ Mq with a star formation rate 
corresponding to a starburst mode can explain basic fea- 
tures of the metal distribution at high z. By accounting 
for the fact that in the presence of UHECR the number 
density of lower massive star-forming halos can be an or- 
der of magnitude higher, and star formation in them sets 
in earlier, even much lower (factor of 20) star formation 
activity than that assumed in Scannapieco et al. (2002) 
can match the observed evolution of metals in the IGM. 
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5. Conclusions 

We have considered the impact of additional ionizing pho- 
tons from UHECR on the lower mass limit of the early 
luminous objects. We showed that in the presence of 
UHECR the early luminous objects are an order of mag- 
nitude less massive, and form earlier in comparison to 
the predictions of the standard recombination history. In 
particular, for e — 0.3 the objects with baryonic mass 
3 X 10^ fif, Mq can form at redshift z = 25. Thus, under 
these conditions the baryons become able to cool and form 
luminous objects at higher redshifts and their fraction in- 
creases by an order of magnitude. A possible consequence 
can be a larger and earlier production of ionizing pho- 
tons in the primary stellar nucleosynthesis required for 
the reionization of the universe. 
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